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Abstract 

From a total data sample of 701.1 pb _1 recorded with e + e~ centre-of-mass energies of 
y^i = 161 — 209 GeV with the OPAL detector at LEP, 11693 W-pair candidate events are 
selected. These data are used to obtain measurements of the W-pair production cross sections 
at 10 different centre-of-mass energies. The ratio of the measured cross sections to the Standard 
Model expectation is found to be: 

data/SM = 1.002 ± 0.011(stat.) ± 0.007(syst.) ± 0.005(theory), 

where the uncertainties are statistical, experimental systematics and theory systematics respec- 
tively. The data are used to determine the W boson branching fractions, which are found to 
be consistent with lepton universality of the charged current interaction. Assuming lepton uni- 
versality, the branching ratio to hadrons is determined to be 67.41 ±0.37(stat.) ±0.23(syst.) %, 
from which the CKM matrix element |V^ S | is determined to be 0.969±0.017(stat.)±0.012(syst.). 
The differential cross section as a function of the W~ production angle is measured for the qqeu 
and qq/iz/ final states. The results described in this paper are consistent with the expectations 
from the Standard Model. 
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1 Introduction 

From 1996 — 2000 the LEP e + e~ collider at CERN operated at centre-of-mass energies, y/s, 
above the threshold for W + W~ production. This paper describes the OPAL measurements of 
the W + W~ production cross section and W branching fractions using this data sample that 
corresponds to a total integrated luminosity of 701.1 pb -1 . The OPAL analysis of W + W~ 
production and decay using data recorded at y/s > 190 GeV has not been published previously. 
For this paper the data recorded at 183 GeV and above have been analysed using the final 
OPAL detector calibration and W pair event selections. The results presented here supersede 
the previous OPAL analysis of the data recorded at y/s = 183 GeV [1] and y/s = 189 GeV [2]. 
The data collected close to the W pair production threshold (y/s = 161 GeV and 172 GeV) have 
not been reanalysed and the corresponding results are described in [3,4]. Furthermore, for the 
reasons explained in Section 3.1, the 183 GeV W + W~— > lv£v data have not been reanalysed 
and the corresponding results are given in [1]. 

In this paper, W + W~ production is defined in terms of the CC03 class [5] of production 
diagrams. These diagrams, which correspond to t-channel v e exchange and s-channel Z/7 
exchange, provide a natural definition of resonant W-pair production. The contributions to the 
event rate from non-CC03 diagrams which lead to the same final states as W-pair production 
(including interference with the CC03 set of diagrams) are treated as additive background. In 
the Standard Model (SM), W + W~ events are expected to decay into fully leptonic (£v£v), semi- 
leptonic (qq^i/), or fully hadronic (qqqq) final states with predicted SM branching fractions of 
10.6%, 43.9% and 45.6% respectively [5]. Here qq denotes a quark and an anti-quark and £v 
denotes a lepton/anti-lepton {£ = e, r) and an anti- neutrino/neutrino. Three separate event 
selections, described in Section 3, are used to identify candidate W + W~ events by their final 
state topologies with £v£v and qq£v candidates classified according to the charged lepton type. 
From the observed event rates in these ten channels (6 £v£v, 3 qqiv and qqqq) measurements 
of the W boson branching fractions and total W + W~ production cross section are obtained. 
The measured branching fraction to hadrons is used to provide a determination of the CKM 
matrix element \V CS \. For the qqev e and qqpv^ decay channels the charge of the W bosons can 
be identified from the charge of the observed lepton. These events are used to determine the 
differential cross section in terms of the W~ polar angle. 
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2 Detector, Data and Monte Carlo 



2.1 The OPAL Detector 

The inner part of the OPAL detector consisted of a 3.7 m diameter tracking volume within a 
0.435 T axial magnetic field. The tracking detectors included a silicon micro-vertex detector, a 
high precision gas vertex detector and a large volume gas jet chamber. The tracking acceptance 
corresponds to approximately | cos 9\ < 0.95 (for the track quality cuts used in this study), where 
9 is the polar angle with respect to the e~ beam direction. The transverse momentum resolution 
for muon tracks is approximately <j Pt /pt = a/(0.02) 2 + (0.0015px) 2 with pr measured in GeV. 
Lying outside the solenoid, the electromagnetic calorimeter (ECAL) consisting of 11704 lead 
glass blocks had full acceptance in the range |cos#| < 0.98 and a relative energy resolution 
for electrons of approximately a E /E ss 0.18 /VE with E measured in GeV. The magnet return 
yoke was instrumented with streamer tubes which served as the hadronic calorimeter. Muon 
chambers outside the hadronic calorimeter provided muon identification in the range | cos#| < 
0.98. Hermeticity for polar angles down to approximately 24mrad was achieved with forward 
detectors designed for measuring electrons and photons. Additional forward scintillator tiles 
were installed in 1998 in order to extend the coverage for detection of minimum ionising particles 
[6]. These forward scintillator tiles were used to improve the Ivlv analysis for the > 189 GeV 
data samples. A detailed description of the OPAL detector can be found in [7]. 

2.2 Data Sample 

From 1996 onwards the centre-of-mass energy of the LEP collider was increased from 161 GeV 
to 209 GeV in several steps. The total integrated luminosity of the data sample considered 
in this paper, evaluated using small angle Bhabha scattering events observed in the silicon 
tungsten forward calorimeter [8], is 701.1 ±2.1 pb _1 . For the purpose of measuring the W + W~ 
cross section these data are divided into ten ^fs ranges listed in Table 1. These ranges reflect 
the main energy steps as the centre-of-mass energy was increased during LEP operation above 
the W + W~ production threshold. 



Range/ GeV 


(v^)/GeV 


r/pb- 1 


160.0- 


-165.0 


161.30 


9.89 


165.0- 


-180.0 


172.11 


10.36 


180.0- 


-185.0 


182.68 


57.38 


185.0- 


-190.0 


188.63 


183.04 


190.0- 


-194.0 


191.61 


29.33 


194.0- 


-198.0 


195.54 


76.41 


198.0- 


-201.0 


199.54 


76.58 


201.0- 


-202.5 


201.65 


37.68 


202.5- 


-205.5 


204.88 


81.91 


205.5- 


-209.0 


206.56 


138.54 


Total 




701.12 



Table 1: The energy binning used for the W + W~ cross section measurements. The \fs range 
covered by each bin, the mean luminosity-weighted value of \fs and the corresponding integrated 
luminosity, C, are listed. 
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2.3 Monte Carlo 



A number of Monte Carlo (MC) samples, all including a full simulation [9] of the OPAL 
detector, are used to model the signal and background processes. For this paper the main MC 
samples for four-fermion final states consistent with coming from the process e + e~ — > W + W~ 
are generated using the KandY [10] program. KandY includes exact 0(pt) YFS exponentiation 
[11] for the W + W~ production process, with 0(a) electroweak non-leading (NL) corrections 
combined with YFS exponentiated 0(a 3 ) leading logarithm (LL) initial state radiation (ISR). 
Final state radiation (FSR) from leptons is implemented in Photos [12] and radiation from 
the quark induced parton-shower is performed by Jetset [13]. The hadronisation within the 
Jetset model is tuned to OPAL data recorded at the Z resonance [14]. For the studies of 
systematic uncertainties the Jetset hadronisation model is compared with the predictions 
from Herwig [15] and Ariadne [16]. 

The KandY generator is also used to produce event weights such that generated events can 
be reweighted to correspond to the CC03 set of diagrams alone. The difference between the full 
set of four-fermion diagrams and the CC03 diagrams alone is used to obtain the four-fermion 
background which includes the effects of interference with the CC03 diagrams. 

The KoralW program [17] is used to simulate the background from four-fermion final states 
which are incompatible with coming from the decays of two W-bosons (e.g. e + e~ — ► qq/i + yU~). 
The two-fermion background processes e + e~ — > Z/7 — > e + e~ — > Z/7 — > t + t~ and 

e + e~ — > Z/7 — > qq are simulated using KK2f [18]. The two fermion process e + e~ — > Z/7 — > 
e + e~ is simulated using Bhwide [19]. Backgrounds from two-photon interactions are evaluated 
using Pythia [20], Herwig, Phojet [21], Bdk [22] and the Vermaseren program [23]. 

The SM predictions for the CC03 e + e~ — > W + W" cross sections above the W + W~ 
threshold region are obtained from the YfsWW [24] and the RacoonWW [25] programs. 
RacoonWW is a complete 0(a) e + e~ — * 4/7 calculation in the double pole approxima- 
tion with ISR treated using a structure function approach. The YfsWW program provides 
the W + W" calculations in KandY. YfsWW and RacoonWW yield nearly identical predic- 
tions for the W + W~ cross sections with an estimated theoretical uncertainty of approximately 
0.5 % [26]. For W-pair production near threshold (the 161 GeV and 172 GeV data) the leading- 
and double-pole approximations used in YfsWW and RacoonWW respectively are no longer 
valid and the predictions are obtained from both calculations using the Improved Born Ap- 
proximation where the theoretical uncertainty is approximately 2%. 

3 e+e- — ► W+W" Event Selection 

The selection of W + W~ events proceeds in three stages, corresponding to the three W + W~ 
decay topologies: W + W~— > Ivlv, W + W~— > qqtv and W + W~^ qqqq- The selections are 
mutually exclusive with only events failing the W + W~— > ivlv selection being considered in 
the W + W~— > qqtv selection, and only events which are not selected as tvtv or qqtv being 
considered for the W + W - ^ qqqq selection. The event selections are essentially unchanged from 
those described in detail in [2] (and references therein) although the W + W~^ iviv selection 
now incorporates features used in the OPAL analysis of di-lepton events with significant missing 
transverse momentum [27]. 

In the centre-of-mass energy range = 161 — 209 GeV, the luminosity-weighted average 
CC03 W-pair selection efficiencies for the £v£v, qqiu and qqqq decay channels are 84%, 84% 
and 86 % respectively. This corresponds to a total efficiency of 85 %. The selection efficiencies, 



5 



broken down into the different lepton ffavours are summarised in Table 2. For the data samples 
away from the W-pair threshold the selection efficiencies depend only weakly on centre-of- 
mass energy. The main features of the selections and associated systematic uncertainties are 
described below in Sections 3.1—3.3. 



bvent 








Efficiencies [%] for W+W" - 








Selection 


evev 


\XV\1V 


TVTV 


evjiv 


evTv 


fiVTV 


qqez/ 


qqfiv 


qqTV 


qqqq 


evev 


7/1 1 
(4.1 


n n 
U.U 


U.o 


0.4 


6.6 


0.1 


0.0 


n n 
U.U 


n n 
U.U 


n n 
U.U 


\xv\xv 


0.0 


77.9 


0.7 


1.4 


0.1 


6.7 


0.0 


0.0 


0.0 


0.0 


TVTV 


0.7 


0.7 


48.1 


0.7 


4.9 


5.6 


0.0 


0.0 


0.0 


0.0 


ev/iv 


2.6 


0.4 


1.4 


76.5 


6.2 


6.9 


0.0 


0.0 


0.0 


0.0 


evTV 


10.3 


0.0 


11.5 


5.6 


64.2 


1.2 


0.0 


0.0 


0.0 


0.0 


fIVTV 


0.2 


9.5 


8.4 


4.3 


0.8 


61.5 


0.0 


0.0 


0.0 


0.0 


qqev 


0.0 


0.0 


0.0 


0.0 


0.2 


0.0 


84.3 


0.1 


4.0 


0.0 


qqpv 


0.0 


0.0 


0.0 


0.0 


0.0 


0.1 


0.2 


88.3 


4.4 


0.1 


qqru 


0.0 


0.0 


0.2 


0.0 


0.0 


0.0 


4.3 


4.4 


61.5 


0.5 


qqqq 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.1 


0.8 


85.9 



Table 2: The luminosity- weighted average selection efficiencies for the CC03 processes for 
y/s = 161 — 209 GeV. The efficiencies include corrections for detector occupancy and tracking 
inefficiencies as described in the text. 



3.1 Selection of W+W~^ lulu events 

The W + W~— > £v£v process results in an event with two charged leptons, not necessarily of 
the same flavour, and significant missing momentum. This characteristic event topology is of 
interest both for measuring aspects of W physics and for exploring the potential production 
of new particles leading to the same experimental signature. The W + W~-^ iviv event selec- 
tion described here first requires events to be selected by the general event selection used by 
OPAL to search for new particles such as pair production of super-symmetric particles which 
decay leptonically [27]. This selection identifies events consistent with there being two charged 
leptons and significant missing transverse momentum. From this sample cuts are applied to 
identify events consistent with being from the W + W~-^ iviv process. This event selection 
takes advantage of changes to the OPAL detector made in 1998. Consequently the data from 
centre-of-mass energies of 161 [3], 172 [4] and 183 GeV [1] have not been reanalysed. 

The general iviv event selection is described in detail in [27] and references therein. The 
selection is formed by requiring that an event be selected by either of two independent event 
selections, referred to in [27] as Selection I and Selection II. Both event selections require ev- 
idence for significant missing transverse momentum and are designed to minimise background 
contributions from SM processes which can lead to an experimental signature of two charged 
leptons and significant missing transverse momentum. In the case of background processes, sig- 
nificant missing transverse momentum can arise from a number of sources: secondary neutrinos 
in tau decays; mis-measurement of the lepton energies and directions; or where high transverse 
momentum particles are incident on poorly instrumented regions of the detector. 

Selection I is designed to retain efficiency for events with low visible energy. Selection II 
is designed for measuring W + W~— > iviv events which usually have substantial visible energy; 
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the selection criteria have been optimised to maximise the statistical power (efficiency multi- 
plied by purity) treating CC03 W + W~— * iviv as signal and SM processes other than iviv as 
background. For both Selection I and Selection II particular care is taken to reject events with 
fake missing momentum due to detector effects. Neither selection attempts to reduce the sensi- 
tivity to non-CC03 sources of iviv events with two detected leptons. There is a large overlap 
in the expected acceptance of the two selections: from the selected MC event sample, 6 % of 
events are selected exclusively by Selection I and 6 % exclusively by Selection II. Conversely, of 
the MC SM background events from processes other than Ivlv, 9 % pass both selections, 32 % 
exclusively pass Selection I and 59 % exclusively pass Selection II. 

Both selections are cut-based and rather involved [27], and only an outline of the main 
points is given here. The most significant variables used are: x min (x max ), the momentum of 
the lower (higher) momentum charged lepton candidate scaled to the beam energy; xt, the 
magnitude of the missing momentum scaled to the beam energy; acop , the supplement of the 
azimuthal opening angle; #™ 1SS , the polar angle of the missing momentum vector; p™ lss , the 
magnitude of the z component of the missing momentum; a miss , the component of the missing 
transverse momentum that is perpendicular to the event thrust axis in the transverse plane; 
and 9f ss = tan- 1 ^" 88 /^ 88 ]. 

Selection I is based on three main requirements: 

• evidence that a pair of charged leptons is produced, where at least one must have pp 
exceeding 1.5 GeV and must satisfy requirements on lepton identification and isolation; 

• evidence of statistically significant missing transverse momentum. For large acoplanarity 
events, acop > tt/2, £t is required to exceed 0.045. For acop < 7r/2, i.e. events where 
the leptons are more back-to-back, a combination of cuts on xt, a™ 188 and #™ 188 is used. 
The cuts depend on the di-lepton identification information; 

• a veto on events with fake missing transverse momentum using the detectors in the forward 
region of the detector. 

Selection I is designed as a general selection for di-lepton events with missing transverse 
momentum. In order to isolate events consistent with the process W + W~— > iviv, additional 
cuts are applied in this analysis to remove events which have relatively low missing transverse 
momentum (an important region for SUSY and other new particle searches but not for W-pair 
production): 

• events are rejected if x max < 0.1; 

• if xt < 0.2, | cos#™ 188 | > 0.7 and x min < 0.3, events are rejected if either x max < 0.15 or 
acop < tt/2 and 6f ss < 0.1; 

• for events with only one reconstructed isolated charged lepton candidate, events are re- 
jected if the net momentum of the additional tracks and clusters not associated to the 
lepton divided by their invariant mass is less than 4. 

Selection II starts from a preselected sample of low multiplicity events and makes little 
use of lepton identification information in the event selection procedure. The first stage of 
the selection is to apply a cone jet-finding algorithm [28] using a cone half-opening angle of 
20° and a jet energy threshold of 2.5 GeV. The majority (90%) of W + W"^ iviv events are 
reconstructed in the di-jet category. For events reconstructed as two jet events, the three most 
important selection criteria are: 



7 



• evidence for missing transverse momentum defined by requiring that xt should exceed 
0.05 by a statistically significant margin; 

• for low acoplanarity events a™ lss should exceed 0.020, primarily to reject events where the 
missing momentum arises from secondary neutrinos from tau decays; 

• a veto on activity in the forward region similar to Selection I. 

Additional selections targeted at three-jet events (often W + W - ^ iviv-f) and single jet events 
(one observed lepton plus evidence for the presence of another lepton) are used to improve the 
overall selection efficiency. 

Events are classified as one of the six possible di-lepton types. For events selected by 
Selection II, the event classification uses both particle identification information and kinematic 
information as described in reference [2]. For events selected exclusively by Selection I the 
di-lepton classification is based on the lepton identification information only. 

3.1.1 W + W _ — > iviv Selection Systematic Uncertainties 

Efficiency Uncertainties: The OPAL trigger and pretrigger systems provide a highly re- 
dundant and efficient trigger for W + W~^ iviv; studies indicate that the trigger inefficiency 
for events selected by these event selections is negligible. The W + W"-^ iviv event selection 
efficiencies are limited mainly by the geometrical acceptance of the detector and the defined 
kinematic acceptance. The latter is implicit in the requirement that the observed final state 
particles have a net visible transverse momentum which significantly exceeds that which could 
be explained by undetected particles at low polar angles. The detector acceptance is well un- 
derstood and factors affecting the kinematic acceptance such as momentum and energy scales 
and resolutions are adequately modelled by the MC simulation. Extensive studies have been 
carried out comparing distributions of the event selection variables in data with MC. In general, 
reasonable agreement is found and quantitative estimates of the individual systematic effects 
are small compared to the statistical errors. In particular, the critical distributions associated 
with requiring missing transverse momentum, such as the a™ ss and the Xt distributions are 
well modelled. As an example, the single most important cut in the two "jet" part of Selection 
II is the cut on a™ lss which leads to a relative loss in the W + W~— * iviv efficiency of 1.1%. 
A conservative estimate of the systematic error on the a^ ms scale of 1% leads to a systematic 
uncertainty of 0.04% on the overall efficiency. As a result of such studies, an overall global 
event selection efficiency systematic uncertainty corresponding to 5 % of the inefficiency prior 
to occupancy corrections is assessed. This systematic uncertainty is taken to be fully correlated 
among centre-of-mass energies and ranges from 0.7% at 189 GeV to 0.8% at 207 GeV. 
Detector Occupancy: The W + W~ — > iviv event selection is sensitive to hits in the various 
sub-detectors which do not arise from the primary e + e~ interaction, termed "detector occu- 
pancy". Backgrounds from the accelerator, cosmic-ray muons, or electronic noise can lead 
to additional hits, energy deposition and even reconstructed tracks being superimposed on 
triggered data events. These detector occupancy effects are simulated by adding to the re- 
constructed MC events the hits, energy depositions and additional "jets" found in randomly 
triggered [29] beam-crossing data events corresponding to the same centre-of-mass energy. The 
detector occupancy corrections are included in the quoted efficiencies of Table 2. They reduce 
the overall efficiency and range from —0.4 % at 189 GeV to —1.0 % at 207 GeV. The variation is 
due to higher beam-related backgrounds at the highest energies. In order to take into account 
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residual deficiencies in the implementation of these post event reconstruction corrections, a 
systematic uncertainty amounting to one half of the correction is assigned. 

The overall Ivlv efficiency systematic uncertainties (for all final states combined) range 
from 0.8% to 1.0 % for centre-of-mass energies of 189 — 209 GeV. 

Background Uncertainties: There are three main sources of background in the W + W~— > Ivlv 
selection: 

• Non-Mi^ Background: Events from processes with no primary neutrinos which man- 
age to fake the missing transverse momentum signature. Important sub- components are 
di-lepton production, in particular tau-pairs, multi-peripheral two-photon processes and 
the four-fermion e + e~ff processes. 

• Non-interfering four-fermion background: Ivlv final states arising from processes 
such as ZZ with primary neutrinos in the final state and with lepton and neutrino flavours 
incompatible with WW production (e.g. fi + fi~v T v T ). 

• Interfering four-fermion background: The Ivlv final states relevant to W + W~— ► Ivlv 
also have significant contributions from diagrams beyond those of CC03 W-pair produc- 
tion, such as Wez/ e , Ze + e~, ZZ and Zv e v e . These contributions, which can also interfere 
with the CC03 diagrams, are treated as an additive background. 

For the centre-of-mass energy range = 161 — 209 GeV, the luminosity-weighted average 
expected background cross sections are listed in Table 3. 



Source of 
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fIVfiV 
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evfj,v 
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qqev 
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17. 


18. 


21. 
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0. 


0. 


0. 
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0. 


0. 
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0. 


0. 


0. 


61. 


3. 
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0. 
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0. 


0. 


0. 


0. 
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0. 


1. 
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0. 
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0. 


0. 
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0. 


0. 


0. 


0. 


1. 


1. 


36. 


0. 
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2. 
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3. 


2. 


1. 
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0. 


qq 


0. 


0. 


0. 


0. 


0. 


0. 


41. 


23. 


78. 


1340. 


e+e-X 


0. 


0. 


7. 


0. 


2. 


1. 


7. 


2. 


3. 


0. 


Total 


23. 


21. 


35. 


23. 


41. 


23. 


152. 


63. 


280. 


1882. 


error 


2. 


3. 


4. 


2. 


3. 


3. 


10. 


5. 


32. 


100. 



Table 3: Luminosity-weighted average background cross sections [fb] in the different event se- 
lection categories. The background cross sections for the qqrv T selection include the corrections 
described in the text. The quoted errors include both statistical and systematic uncertainties. 

The overall systematic uncertainties on the background cross sections for each di-lepton 
class and at each centre-of-mass energy are calculated by summing up the contributions in the 
following categories. The uncertainties within each category are assumed to be fully correlated 
among di-lepton channels and centre-of-mass energies. 
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• For events from di-lepton production the theoretical uncertainties are negligible. In this 
case it is simulation of the detector response that dominates the uncertainty on the 
background. Events are selected due to either mis-measurements of the variables used in 
the selection or from the tails of the t + t~ decay distributions. An overall background 
systematic uncertainty of 10 % is assessed. 

• A 5 % systematic uncertainty is assigned to the background expectations from genuine 
iviv events coming both from non-interfering four-fermion background final states and 
from the non-CC03 contribution to final states where the four fermions are compatible 
with being from W-pair production. 

• A 10 % systematic uncertainty is assigned to the background expectations from e + e~ff 
and the remaining small contributions from other four-fermion processes, reflecting the 
theoretical error on simulation of processes like Ze + e~. 

• For events from the multi-peripheral e + e~ — > e + e~X process an uncertainty of 30 % is 
assigned. The uncertainty reflects the size of the discrepancy in the modelled number 
of events exclusively rejected using the forward scintillating tiles, a category of events 
dominated by multi-peripheral backgrounds. 

Event Classification Uncertainties: There are two aspects to the di-lepton flavour classifi- 
cation of selected W + W~^ iviv candidates. Firstly, the algorithms for leptons to be identified 
as electrons, muons or hadronically decaying taus. These make use of many of the techniques of 
lepton identification used by OPAL in studies at the Z. Secondly, the kinematic re-classification 
algorithm based on scaled momentum which re-classifies soft leptons identified as electrons or 
muons as probable secondary leptons from taus, and uses electromagnetic calorimeter and muon 
information to re-assess whether highly energetic leptons initially not identified as electrons or 
muons are more consistent kinematically with prompt electrons or muons. The classification 
efficiency systematic uncertainty for genuine electrons and muons is assessed to be 2 % based 
on the understanding of the lepton identification information in the large e + e~ — > i + i~ samples 
recorded at LEP1. The kinematic re-classification, which relies mainly on measurement of the 
lepton energy, reduces the systematic uncertainties on the efficiencies for the individual final 
state lepton channels to the 1 % level. In the extraction of the SM parameters that follows it has 
been verified that the effects of the iviv classification systematic uncertainties are small. Nev- 
ertheless, the effects of the classification systematic uncertainties and correlations are included 
in the analysis. 

3.1.2 W+W-^ iviv Results 

Using the KandY MC samples the luminosity- weighted average CC03 W + W~^ iviv event 
selection efficiency in the 189— 209 GeV centre-of-mass energy range is estimated to be (84.7 ± 
0.8) %. The inclusive selection efficiencies for the different centre-of-mass energies are listed 
in Table 4. The efficiencies for the different final states depend mostly on the number of 
taus present. The luminosity-weighted average efficiencies are 89.4%, 83.2% and 71.9% for 
final states with zero, one and two taus respectively. For the 189—209 GeV data the selection 
efficiency does not depend strongly on centre-of-mass energy. The luminosity weighted effi- 
ciencies of the W + W~— > iviv selection for the individual channels are given in Table 2. The 
efficiencies/numbers of expected events in all tables include the detector occupancy corrections 
described above. 
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In total, 1188 events are selected as W + W~^ iviv candidates compared to the SM expec- 
tation of 1138 ±9 (the numbers refer to the entire data set from 161 — 209 GeV). Figure 1 shows 
kinematic distributions for reconstructed W + W~— > iviv event samples. The data distributions 
are in good agreement with the MC expectations. The numbers of selected iviv events at each 
energy are used to determine the cross sections for e + e~ — > W + W~ — > iviv given in Table 4. 
The measured cross sections are in agreement with the SM expectations. 
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1.80 


204.88 


81.9 
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83.5± 1.0 


16.3± 0.5 


1.82± 0.17± 0.02 


1.81 


206.56 


138.5 


239 
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1.81 



Table 4: Measured cross sections for the CC03 process e + e~ — > W + W~ — > iviv. For the iviv 
selection the data below ^/s = 188.63 GeV have not been reanalysed and the results are taken 
from [1,3,4]. The errors on the cross sections are statistical and systematic. The numbers of 
selected events, the iviv selection efficiencies and the expected numbers of background events 
are also listed. The backgrounds include a small contribution from semi-leptonic W + W" decays 
which for the cross sections are taken to be fixed to their SM expectations. 



3.2 Selection of W + W — > qqjtv events 

The W + W~— > qqiv selection consists of three separate selections, one for each type of semi- 
leptonic decay. Only those events which are not already selected as iviv candidates are consid- 
ered by these selections. For each of the W + W~— > qqev e , W + W"— > qqpv^, and W + W"-^ qqrv T 
event selections, the main part is a relative likelihood method to reject the potentially large 
e + e~ — > qq background. In the first stage, the W + W"— > qqev e and W + W~-^ qqp^n likelihood 
selections are performed. The W + W~^ qqrv T likelihood selection is only applied to those 
events which have not already been selected. Finally, events passing either the W + W~^ qqez/ e 
or the W + W~— >• qq^v^ selections may then be reclassified as W + W~— > qqrv T candidates. 

The W + W~— > qqiv event selections used here are almost identical to those described in 
previous OPAL publications [1,2]. However, using the entire OPAL W + W~ data has resulted 
in an improved understanding of the selection efficiencies and backgrounds. Using the improved 
estimates of the systematic uncertainties, the cut on the relative likelihood variable used to 
select qqrv T candidates was re-optimised to minimise the total uncertainty (statistical and 
systematic) for this channel. As a result the cut on the likelihood was raised from 0.5 to 0.8 
which reduces the efficiency by about 5%. This loss in efficiency is more than compensated by 
the factor of two reduction in background and the corresponding reduction in the associated 
systematic uncertainties. 
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3.2.1 Event Selection 



The W + W~— > qqiv event selection utilises the distinct topology of W + W"-^ qqiv events; 
missing energy and a high energy (usually isolated) lepton. The selection consists of six stages, 
which can be summarised as: 

• loose preselection: a loose preselection to remove events with low multiplicity or little 
visible energy. 

• lepton candidate identification: identification of the observed track in the event which 
is most consistent with being from the leptonic decay of a W boson. Candidate lepton 
tracks are identified for each of the qqez/ e , qquz/^ and qqrz/ T hypotheses. 

• preselection: different sets of cuts are applied for W + W~ — >■ qqez/ e , W + W~^ qq/i^, 
and W + W - -^ qqTu T to remove events clearly incompatible with being signal (e.g. events 
are rejected if the total visible energy in the event is less than 0.3 of the centre-of-mass 
energy). 

• relative likelihood selection: different relative likelihood selections are used to iden- 
tify W + W~— > qqez/ e , W + W~— > qqu^ M , and W + W - ^ qqTv T candidates. The probability 
density functions used in the likelihood selections are obtained from MC at the different 
centre-of-mass energies. The variables used are either related to the properties of the lep- 
ton candidate (e.g. the lepton energy and degree of isolation) or the kinematic properties 
of the event (e.g. the total visible energy and the magnitude of the missing momentum). 

• decay classification: identification of qqrv T candidates from events which were origi- 
nally selected as qqez/ e or qqjiu^. 

• four-fermion background rejection: rejection of four-fermion backgrounds qq£ + £~, 
Wez/ e , Ze + e~ and qqvu. 

The first four stages, described in detail in [4], are optimised for the rejection of the e + e~ — > qq 
background which, for the centre-of-mass energies considered here, has an expected cross sec- 
tion of between four and seven times larger than the W-pair production cross section. The 
most important feature of the selection is the looseness of the identification of possible lepton 
candidates. For both the W + W"^ qqeu e and W + W~— * qq/ii^ selections the track which is 
most consistent with being from a leptonic W-decay is identified. The lepton track identifica- 
tion is based on an absolute likelihood taking into account momentum, isolation and lepton 
identification variables. To avoid associated systematic uncertainties only very loose cuts are 
placed on the lepton identification likelihood. The lepton identification likelihood is then used 
as one of the input variables in the likelihood event selection. In this way the presence of 
either a good isolated lepton candidate or significant missing transverse momentum is usually 
sufficient for an event to be selected. This redundancy leads to high efficiency and reduces the 
dependence of the selection on the detailed simulation of the events and, consequently, leads 
to relatively small systematic uncertainties. 

Because of the limited use of lepton identification information, approximately 33 % of 
W + W~— > qqrv T events are accepted by at least one of the qqez/ e and qq/xz^ likelihood se- 
lections. In addition, approximately 4% of the W + W~— > qqeu e and W + W~— * qq/xz^ events 
pass both the qqeu e and qq/xz^ likelihood selections. Such events usually result from there 
being a genuine electron from a W-boson decay and a track from one of the jets being tagged 
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as muon-like, or vice versa. Consequently additional likelihood selections, based primarily on 
lepton identification variables and track momentum, are used to categorise events passing the 
qqez/ e and qq/if^ likelihood selections into the three possible leptonic W-decay modes. The 
largest systematic uncertainties in the efficiencies for selecting W + W~— > qqtv events are asso- 
ciated with this step. 

Only events which failed the W + W~ — > qqeu e and W + W~^ qq/^ likelihood are passed to 
the W + W - ^ qqT^r event selection. The W + W - ^ qqTu T event selection consists of separate 
selections for four possible tau decay signatures: r — > euu, t — > fiuu, single prong hadronic 
decay modes and three prong hadronic decay modes. The main difference between these selec- 
tions is the power of the variables used to identify possible tau decay products and the relative 
level of backgrounds. An event is considered a qqru T candidate if it passes any one of these 
four selections. 

Because the W + W~— > qqlv likelihood selections are designed to reject the dominant e + e~ — > 
qq background they have a significant efficiency for other four-fermion processes, e.g. qqez/ e 
final states produced by the single W (Wez/ e ) diagrams and qq£ + £~ production (mainly via 
e + e~ — > ZZ). Additional four-fermion background rejection cuts are applied to events pass- 
ing the likelihood selections to reduce backgrounds from these processes. The four-fermion 
background rejection consists of three separate parts. Cuts are applied to selected qqez/ e and 
qq/iz/ M candidates to reduce backgrounds from qqe + e~ and qq/i + /i~ final states where both 
leptons are observed in the detector. Because of the lack of a clear signature for a lepton 
in W + W~^ qqTU T events, the selection places more weight on missing transverse energy to 
reject e + e~ — > qq. Consequently the W + W~^ qqru T selection accepts approximately 40% of 
hadronically decaying single W events (Wez/ e — > qqeu e ). In these events the electron is usually 
produced in the far forward region beyond the experimental acceptance and a fragmentation 
track is mis- identified as a r lepton decay product. To reduce this background, an additional 
likelihood selection is applied which separates W + W~— >■ qqru T from Wez/ e . This also rejects 
background from e + e~ — > qquV. Background in the W + W"^ qqez/ e selection from the Ze + e~ 
final state, where the Z decays hadronically and one electron is far forward, is reduced with two 
kinematic fits, the first using the hypothesis that the event is W + W~— > qqeu e and the second 
using the Ze + e~ hypothesis. 

In addition to the likelihood selections, cut based selections are used to identify W + W~-^ qqez/ e 
and W + W~— > qq//i^ events where the lepton track is either poorly reconstructed or is beyond 
the tracking acceptance. These 'trackless' selections require clear evidence of an electron or 
muon in the calorimeter or muon chambers consistent with the kinematics of a W + W~^ qqiu 
event, without explicitly demanding a reconstructed track. These additional selections improve 
the overall efficiency by approximately 3 % (5 %) for W + W~^ qqei/ e (W + W~— > qq/iz^) events, 
and more importantly result in a reduction in the systematic uncertainties associated with the 
modelling of the forward tracking acceptance. 

3.2.2 Systematic uncertainties 

Table 5 lists the various contributions to the systematic uncertainty on the qqez/ e , qq/xz^ and 
qqrz/ T selection efficiencies. Many of the potential systematic effects primarily affect the clas- 
sification of selected qq£z/ events rather than the overall qq£z/ efficiency. Amongst the effects 
studied were: 

i) Finite MC statistics of the KandY MC samples used to determine the efficiencies. 

ii) The fragmentation and hadronisation systematic uncertainties are studied with fully simu- 
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Signal efficiency error (%) 
Event Selection W+W~ -) 



Source of uncertainty 


qqez/ e 


qq/iz/ M 


qqrv T 


qqiv 


i) 


MC Statistics 


0.07 


0.06 


0.10 


0.04 


ii) 


WW Fragmentation 


0.25 


0.20 


0.50 


0.20 


iii) 


Tau candidate ID 


— 


— 


0.60 


0.20 


iv) 


0(a) QED/Electroweak 


0.09 


0.05 


0.03 


0.04 


v) 


ISR and FSR 


0.07 


0.12 


0.10 


0.03 


vi) 


ECAL energy response 


0.11 


— 


0.08 


0.03 


vii) 


Track momentum response 


0.07 


0.05 


0.08 


0.02 


viii) 


Jet energy response 


0.01 




0.02 


0.01 


ix) 


Tracking Losses 


0.30 


0.05 


0.06 


0.10 


x) 


Detector Occupancy 


0.03 


0.03 


0.06 


0.03 


xi) 


Preselection 


0.10 


0.10 


0.15 


0.12 


xii) 


Likelihood Selection 


0.30 


0.10 


0.40 


0.10 




Other 


0.04 


0.03 


0.02 


0.03 


Total 


0.54 


0.30 


0.91 


0.36 



Table 5: Sources of uncertainty on the W + W~— > qqiv selection efficiencies. The errors quoted 
apply to the selection efficiency for the combined yfs = 183 — 209 GeV data set. Entries where 
the systematic error estimate is less than 0.01 % are denoted by — . The errors on the combined 
qqiv selection take into account correlations between the separate channels. 

lated MC W + W~^ qqiv samples where the hadronisation process is modelled using Jetset, 
Herwig or Ariadne. In addition, the parameters a q , b, Aq C d, and Q Q of the Jetset frag- 
mentation model are varied by one standard deviation about their tuned values [14]. 

iii) The largest single systematic uncertainty in the qqiv selection is due to an identified defi- 
ciency in the MC simulation of isolated tracks from the fragmentation/hadronisation process. 
Such tracks, if sufficiently isolated can have similar properties to those from hadronic tau de- 
cays. In data there is a clear excess of low momentum tracks which have been identified as 
the best tau decay candidate compared to the MC expectation. This excess persists at all 
stages in the event selection; for example, there is a ~10 % excess of data events passing the 
W + W - -^ qqrv T preselection cuts (a sample dominated by background from e + e~ — > qq). To 
assess the impact on the qqrv T analysis, a control sample of two jet events is formed by remov- 
ing the tracks and calorimeter clusters associated with the lepton in selected qqev e and qq^v^ 
events. The full qqrv T event selection is applied to these events and the selection efficiency is 
found to be 7.3 ± 4.6 % higher in data than the MC expectation. Again there is a clear excess 
(25 ±7%) of isolated tracks with momenta less than 5 GeV. This data sample is used to pro- 
vide a momentum dependent correction factor which is used to reweight all MC events where a 
fragmentation track is identified as the best tau candidate. After applying this correction, the 
data/MC agreement at all stages in the qqrv T selection is significantly improved. The effect of 
this correction is to increase the expected background from qqvv and single-W (Wev e ) events. 
Because qqrv T events can also be selected on the basis of a fragmentation track, the predicted 
selection efficiency for qqrv T events is also increased by 0.6 %. The full size of the corrections to 
efficiency and background are assigned as (correlated) systematic errors in the qqrv T selection. 
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iv) The selection efficiencies are sensitive to hard photon radiation in the W-pair production 
process. The OPAL data are consistent with the predictions from KandY [30]. Potential 
systematic biases are estimated by reweighting the KandY MC samples so as to turn off the 
0(a) electroweak treatment of radiation from the W-bosons. 

v) A conservative estimate of the possible biases arising from FSR from the lepton or tau decay 
products is investigated by reweighting the MC so as to change the rate of such FSR by ±50 %. 
This mainly affects the classification of selected events. The selection efficiencies are found to 
be insensitive to the detailed treatment of ISR. 

vi) , vii) and viii) Uncertainties in the detector calibration, linearity of energy response and 
MC simulation of the energy resolution were studied in detail for the OPAL analysis of the W- 
boson mass [31]. The uncertainties related to ECAL energy, track momentum and jet energy 
response described therein are propagated to the event selection. 

ix) Z — > events are used to study the tracking efficiency for electrons and muons. It is 
found that the MC overestimates the efficiency for reconstructing electron and muon tracks in 
the forward region, | cos#| > 0.9. The effect on the selection efficiency is reduced by a factor of 
approximately three due to the trackless selections. The MC efficiency estimates are corrected 
and the full size of the correction is assigned as a systematic error. 

x) Randomly triggered events recorded throughout the data-taking period are used to assess the 
impact of energy deposits in the detector (particularly in the forward luminosity calorimeters) 
which can result in the event being vetoed. As a result, the MC efficiencies were corrected and 
half the correction assigned as a systematic uncertainty. 

xi) The event preselection cuts remove approximately 1 % of qqtv events. Possible systematic 
effects specifically associated with the preselection (in addition to those described above) are 
studied applying the likelihood selection to all events failing just one of the preselection cuts. 
There is no evidence of any systematic bias and the statistical precision of the study is used to 
assign the systematic uncertainty. 

xii) The MC expectation for each of the variables used in the likelihood selection is compared 
to the observed distribution for the selected events. The ratio of data to MC is used to define 
bin-by-bin corrections for each distribution. These corrections are propagated back into the 
likelihood selection and the associated systematic errors are obtained from the resulting changes 
in the selection efficiencies. 

Background Uncertainties: Table 3 shows the background cross sections and total uncer- 
tainties for the three qqiv selections. The largest contributions to the background in the qqlv 
selections are from the four fermion final states qqez/ e , qq£ + £~ and qqvV and from e + e~ — > qq. 
In the qqrv T selection, the uncertainties on the four fermion backgrounds are dominated by the 
correction for isolated low momentum tracks described above. The qqez/ e background mainly 
arises from the single W process (including interference with the CC03 diagrams); a 5% un- 
certainty on this cross section is assumed [26]. Background from the e + e~ — > qq process mainly 
arises from radiative return events with an unobserved photon in the beam direction where 
a hadronisation track is mis-identified as the lepton. The e + e~ — > qq background is assigned 
a 10 % systematic uncertainty for the MC modelling of the hadronisation process (based on 
comparisons of Pythia, Herwig and Ariadne). The MC estimate of this background rate is 
checked using control samples constructed from the data directly. For the background, 'fake' 
events are constructed by boosting hadronic Z events recorded at \/s = 91 GeV to the invari- 
ant mass distribution expected of quark pairs at the appropriate y/s. There is an additional 
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11% uncertainty on the e + e~ — > qq background in the qqez/ e selection from uncertainties in 
the rate at which high energy photon conversions fake an electron. The backgrounds from 
multi-peripheral two photon processes (almost entirely from hadronic final states rather than 
from e + e~ — > e + e~l + l~) are assigned a systematic uncertainty of 50% to cover the variation 
in predictions obtained from different generators. 

3.2.3 W+W-^ qq£v Results 

Using the KandY MC samples the inclusive qqlv selection is estimated to be 83.8 ± 0.4% 
efficient for W + W - ^ qqlv events. The selection efficiencies for the different centre-of-mass 
energies are listed in Table 6. Above the W + W~ threshold region the selection efficiency 
does not depend strongly on the centre-of-mass energy. The luminosity weighted efficiencies 
of the W + W~-^ qqlv selection for the individual channels are given in Table 2. The effi- 
ciencies/numbers of expected events in all tables include small corrections (0.1 — 0.3%) which 
account for tracking losses which are not modelled by the MC simulation of the OPAL detector. 
The effect of detector occupancy from beam-related backgrounds is also included as is the small 
correction associated with the identification of tau candidates described above. 

In total 4572 events are selected as inclusive W + W~-^ qqlv candidates in agreement with 
the SM expectation of 4622 ± 28. Figure 2 shows distributions of the reconstructed energy of 
the lepton in the qqev e , qq/iv^, and qqrv T selections and the summed distribution. The data 
distributions are in good agreement with the MC expectations. 

The numbers of selected qqlv events at each energy are used to determine the cross sections 
for e + e~ — > W + W~ — > qqlv given in Table 6. The results are obtained assuming the small 
backgrounds from Ivlv and qqqq are given by the SM. The measured cross sections are in 
agreement with the SM expectations. 





C 


N 


Efficiency 


Background 


a(W + W- -> qqlv) 


SM 


[GeV] 




[events] 


[%] 


[events] 


[pb] 


[pb] 


161.30 


9.9 


12 


63.6± 2.5 


1.4± 0.5 


1.68± 0.55± 0.07 


1.58 


172.11 


10.4 


55 


84.2± 1.0 


4.6± 0.8 


5.77± 0.85± 0.07 


5.31 


182.68 


57.4 


357 


84.2± 0.4 


22. 1± 2.1 


6.93± 0.39± 0.05 


6.74 


188.63 


183.0 


1171 


84.6± 0.4 


89.8± 5.7 


6.98± 0.22± 0.05 


7.13 


191.61 


29.3 


176 


84.6± 0.4 


15. 1± 1.0 


6.48± 0.54± 0.05 


7.26 


195.54 


76.4 


554 


84.1± 0.4 


43.6± 2.6 


7.94± 0.37± 0.05 


7.38 


199.54 


76.6 


494 


83. 7± 0.4 


44.8± 2.7 


7.01± 0.35± 0.05 


7.46 


201.65 


37.7 


255 


83.6± 0.4 


22. 1± 1.3 


7.39± 0.51± 0.05 


7.48 


204.88 


81.9 


523 


83.9± 0.4 


52.3± 3.2 


6.85± 0.33± 0.05 


7.50 


206.56 


138.5 


975 


83.6± 0.4 


86.9± 5.1 


7.67± 0.27± 0.05 


7.51 



Table 6: Measured cross sections for the process e + e~ — > W + W~ — * qqtv. For the qqlv selection 
the data below y/s = 182.68 GeV have not been reanalysed and the results are taken from [3,4]. 
The errors on the cross sections are statistical and systematic respectively. The numbers of 
selected events, qqlv selection efficiencies and expected numbers of background events are also 
listed. The backgrounds include fully-leptonic and fully-hadronic W + W" decays for which the 
cross sections are taken to be their SM expectations. 



16 



3.3 Selection of W + W — > qqqq events 

The selection of fully hadronic W + W~— > qqqq events is performed in two stages using a cut- 
based preselection followed by a likelihood selection procedure. This likelihood selection is 
primarily designed to reject the dominant background from the e + e~ — > qq process where 
the di-quark system fragments into a four jet topology. No attempt is made to discriminate 
against the neutral current process ZZ — > qqqq for which the cross section is at least an order 
of magnitude smaller than that for W + W~^ qqqq- The preselection and likelihood selection 
variables are unchanged from those described in previous OPAL publications [2] although the 
tuning of the likelihood discriminant is updated for different ranges of y/s. 

3.3.1 Event Selection 

All events which are classified as hadronic [32] and which have not been selected by either the 
£i>£i> or the qqiu selections are considered as candidates for the W + W~— > qqqq selection. In 
addition, any event which is identified and rejected as a four-fermion background event in the 
qqiv selection is also rejected as a qqqq candidate event. 

Tracks and calorimeter clusters are combined into four jets using the Durham algorithm [33] 
and the total momentum and energy of each jet is corrected for double-counting of energy [34]. 
To remove events which are clearly inconsistent with a fully hadronic W + W~ decay, candidate 
events are required to satisfy a set of preselection cuts including a cut on minimum visible 
energy (70% of y/s), minimum invariant mass (75% of y/s), and minimum multiplicity per jet 
(one track). The most important preselection cut is log 10 (W-42o) < [35], where W420 is the 
QCD matrix element calculated as an event weight formed from the tree level 0(a 2 s ) matrix 
element [36] for the four jet production processes (e + e~ — > qq — > qqqq, qqgg). The value of 
W420 is determined by using the observed momenta of the four reconstructed jets as estimates 
of the underlying parton momenta which are input to the matrix element calculation. The best 
discriminating power between signal and background was found using a variable defined as the 
largest value of the W420 matrix element from any of the 24 possible jet-parton associations in 
each event. 

The preselection requirements reject around 95 % of the e + e~ — > qq events which comprise 
the dominant source of background in the W + W~^ qqqq event selection, while the preselection 
efficiency for the hadronic W + W~^ qqqq decays is estimated to be 90 — 93% depending on 

Events satisfying the preselection cuts are classified as signal or background based upon 
a four variable likelihood selection. The following likelihood variables are selected to provide 
a good separation between the hadronic W + W - ^ qqqq signal and the e + e~ — > qq four jet 
background, while minimising the total number of variables used: 

• log 10 (W42o), the QCD four jet matrix element; 

• log 10 (Wcco3), the Excalibur matrix element [37] for the CC03 process (W + W~^> qqqq); 

• log 10 (y45), the logarithm of the value of the Durham jet resolution parameter at which an 
event is reclassified from four jets to five jets; 

• event sphericity. 

Figure 3 shows the distribution of these four likelihood variables for all preselected events found 
in the 183—209 GeV data. To improve the statistical power of this selection, a multi-dimensional 
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likelihood technique is used to account for the correlations between the four likelihood input 
variables [38]. Most of the separation between the signal and background events is provided 
by the two matrix element values log 10 (W / cco3) an d log 10 (W42o), which is related to the relative 
probability that the kinematics of the observed event are consistent with signal or background 
production respectively. While the likelihood input variables are the same for events in all y/s 
ranges, the likelihood discriminant functions are separately calculated from CC03 signal and 
e + e~ — > qq background MC samples in three ranges of y/s: 185 — 194 GeV, 194 — 202.5 GeV, 
and 202.5 — 209.0 GeV. Candidate events at y/s below 185 GeV are unchanged from previous 
OPAL publications [1,3,4]. 

An event is selected as a hadronic W + W~ — > qqqq candidate if the likelihood discriminant 
variable, also shown in Figure 3, is greater than 0.4. This cut value was chosen to maximise 
the expected statistical power of this selection assuming the SM rate for CC03 production. 

3.3.2 Background Estimation 

The accepted e + e~ — > qq background is estimated from KK2f MC samples, with Pythia 
Herwig and Ariadne hadronisation being used as cross-checks. To reduce the uncertainty on 
this background estimate, a technique to measure this rate directly from the data is used. By 
comparing the number of events seen in data and MC in the range < log 10 (W4 2 o) < 1 which 
would otherwise pass the preselection cuts, the overall four jet background rate predicted by the 
MC is normalised to the observed data. This procedure is performed and applied separately 
in the three y/s selection ranges described above. A luminosity-weighted average correction 
over the full y/s range of (—1.4 ± 1.7)% is found for the default KK2f samples, where the 
uncertainty is the statistical precision of the normalisation procedure. The observed data and 
corrected MC expectation in this sideband background region are shown in Figure 3. The 
expected contamination from CC03 production in this region is less than 3%, resulting in a 
negligible bias on the extracted CC03 cross section. 

3.3.3 Selection Uncertainties 

The main systematic uncertainty on the selection efficiency results from the modelling of the 
QCD hadronisation process. This uncertainty is estimated by comparing the selection efficiency 
predicted using the Jetset hadronisation model with alternative models including Herwig, 
Ariadne and an older version of the OPAL Jetset tuning [39]. These variations cover the 
observed data/MC differences such as the 2/45 distribution shown in Figure 3. The uncertainty in 
the selection efficiency from the modelling of the hadronisation process is almost exclusively due 
to the preselection requirements, and is found to be independent of y/s. The largest observed 
deviation in selection efficiency is taken as the systematic uncertainty, resulting in an estimated 
relative uncertainty of 0.9% which is fully correlated between different y/s samples. 

Cross-checks of this uncertainty are performed by comparing the observed shapes of both 
the preselection and selection variables seen in data to those predicted by the signal MC sam- 
ples. After subtracting the expected background, the differences between observed data and 
expected MC signal distributions are comparable to the variations observed within the different 
hadronisation models themselves. In addition, the effect of directly varying the parameters <r q , 
b, Aqcd, and Qq of the Jetset hadronisation model by one standard deviation about their 
tuned values [14] as was done for previous OPAL results [2] leads to similar uncertainties. 

Additional uncertainties on the modelling of the underlying hard process are evaluated by 
comparing CC03 events produced by KandY with other generators (Excalibur, Pythia, and 
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grc4f [40]). Uncertainties on the detector modelling are evaluated from direct comparison of 
data distributions with MC predictions, and are generally smaller than the observed differences 
seen between the different hadronisation models. Possible biases related to final state interac- 
tions between the hadronic systems produced by different W bosons have been evaluated for 
colour-reconnection effects [41] and Bose-Einstein correlations [42]. These effects are found to 
be small, and the total change in predicted selection efficiency when these effects are included 
in the hadronisation model is taken as the systematic uncertainty. 

3.3.4 Background Uncertainties 

The dominant uncertainty on the expected background rate comes from the modelling of the 
hadronisation process, particularly in e + e~ — > qq events. This uncertainty is evaluated in the 
same manner as the hadronisation uncertainty for the signal efficiency, using large MC samples 
produced with a variety of hadronisation models, and taking the largest observed deviation 
as an estimate of the systematic uncertainty. The background normalisation procedure has 
been consistently applied during these systematic checks. The uncertainty on the estimated 
background is about 75 fb (the exact value depends on the centre-of-mass energy) which is 
taken to be fully correlated between different ^/s samples. The uncertainty from modelling of 
the hadronisation process for the background estimation is found to be largely uncorrelated 
with the uncertainty on the signal efficiency. 

The background normalisation procedure contributes an additional, statistical uncertainty 
to the background estimation of about 3 % which is uncorrelated between different a/s ranges. 
Additional uncertainties in the non-CC03 four-fermion background are estimated by comparing 
the expectations of KoralW, grc4f, and Excalibur. This background is predominantly 
from the neutral current process ZZ — > qqqq, of which only 20% is in final states with direct 
interference with the CC03 diagrams. In each case, the single largest difference observed in a 
set of systematic checks is taken as an estimate of the uncertainty. 

3.3.5 W+W-^ qqqq Results 

The luminosity- weighted efficiency of the likelihood selection for W + W~— > qqqq events is es- 
timated from KandY MC samples to be 85.9 ± 0.9 %, where the error represents an estimate 
of the systematic uncertainties. A total of 5933 W + W"^ qqqq candidate events are selected 
compared to the expectation of 5845.2 ± 67.5. The luminosity-weighted purity of the selected 
event sample is 77%. The selection efficiencies for the different centre-of-mass energies are 
listed in Table 7. For the 189—209 GeV data the selection efficiency does not depend strongly 
of centre-of-mass energy. The numbers of selected qqqq events at each energy are used to deter- 
mine cross sections for e + e~ — > W + W~ — ► qqqq, also listed in Table 7. The results are obtained 
assuming the small backgrounds from iviv and qq£v are given by the SM. The measured cross 
sections are in agreement with the SM expectations. 

4 Measurement of the W + W~ cross section 

The observed numbers of selected W + W~ events are used to measure the W + W~ production 
cross section and the W decay branching fractions to leptons and hadrons. The measured cross 
section corresponds to that of W-pair production from the CC03 diagrams as discussed earlier. 
The expected four-fermion backgrounds quoted throughout this paper include contributions 
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L 


N 


Efficiency 


Background 


^(W+W- -> qqqq) 


SM 


[GeV] 


[pb- 1 ] 


[events] 


[%] 


[events] 


[pb] 


[pb] 


161.30 


9.9 


14 


56. 7± 3.5 


3.4± 0.4 


1.88± 0.67± 0.14 


1.64 


172.11 


10.4 


54 


70.3± 3.0 


13.1± 1.9 


5.62± 1.01± 0.24 


5.52 


182.68 


57.4 


439 


86.3± 0.9 


98.1± 6.8 


6.89± 0.42± 0.11 


7.00 


188.63 


183.0 


1553 


86.6± 0.9 


339.5± 17.8 


7.66± 0.25± 0.12 


7.41 


191.61 


29.3 


245 


86.2± 0.9 


55.2± 2.8 


7.51± 0.62± 0.12 


7.54 


195.54 


76.4 


709 


87.2± 0.9 


152. 6± 7.8 


8.35± 0.40± 0.12 


7.67 


199.54 


76.6 


643 


86. 7± 0.9 


150.6± 7.7 


7.42± 0.38± 0.11 


7.75 


201.65 


37.7 


342 


86.6± 0.9 


75.8± 3.8 


8.16± 0.57± 0.12 


7.77 


204.88 


81.9 


683 


86.3± 0.9 


159.9± 8.2 


7.40± 0.37± 0.11 


7.79 


206.56 


138.5 


1251 


86. 1± 0.9 


274.4± 13.9 


8.19± 0.30± 0.12 


7.80 



Table 7: Measured cross sections for the process e + e~ — > W + W~ — > qqqq. For the qqqq 
selection the data below y/s = 182.68 GeV have not been reanalysed and the results are taken 
from [3,4]. The errors on the cross sections are statistical and systematic respectively. The 
numbers of selected events, qqqq selection efficiencies and expected numbers of background 
events are also listed. The backgrounds include fully-leptonic and semi-leptonic W + W~ decays 
which for the cross sections are taken to be fixed to their SM expectations. 

from both non-CC03 final states and the effects of interference with the CC03 diagrams. 
Mis-identified CC03 final states are not included in the background values listed in Table 3, 
but rather are taken into account by off-diagonal entries in the efficiency matrix. Table 8 
summarises the event selections in the ten W + W~ decay topologies. 

The W + W~ cross section and branching fractions are measured using data from the ten 
separate decay channels. The physical parameters (cross sections, branching ratios, etc.) are 
obtained from fits where all correlated systematic uncertainties are taken into account. The 
total cross section is obtained from a maximum likelihood fit to the numbers of events in 
the ten decay channels from data at all centre-of-mass energies allowing the cross sections 
at each centre-of-mass energy to vary and assuming the SM branching fractions. Efficiency, 
background, and luminosity systematic uncertainties are included as nuisance parameters with 
Gaussian penalty terms in the likelihood function [43]. Correlations are accounted for in the 
covariance matrix of the nuisance parameters associated with the systematic uncertainties. The 
results are listed in Table 9 and shown in Figure 4. In both cases the results are compared to 
the SM expectation which is taken to be the mean of the cross sections predicted by YfsWW 
and RacoonWW (on average the predicted cross section from YfsWW is 0.2 % higher than 
that from RacoonWW). The results do not differ significantly if the SM branching fractions 
are left unconstrained in the fit. When compared to the SM expectations, the 10 cross section 
measurements in Figure 4 yield a x 2 of 15.5 (11 % probability). When the 100 individual event 
counts used to obtain the cross sections (ten channels x ten y^s bins) are compared to the SM 
expectation the x 2 obtained is 94.5 for 100 degrees of freedom. The OPAL W + W~ data are 
consistent with the SM expectation. The cross sections listed in Table 9 differ from than the 
sums of the exclusive cross sections from the separate channels (listed in Tables 4, 6 and 7) 
because of the constraint to the SM branching ratios and the larger systematic errors and in 
the qqqq channel. 

A fit to the data where the expected cross sections at all centre-of-mass energies are given 
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Selection 
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Data/Expected 
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TVTV 


71.8% 


79.5% 


122. 2± 3.4 


131 


1.072+0.094+0.028 


evfj,v 


91.8% 


93.9% 


264.8± 3.2 


251 


0.948+0.060+0.012 


evrv 


81.9% 


88.5% 


250.5± 4.2 


256 


1.022+0.064+0.017 


fIVTV 


75.6% 


92.6% 


220.9± 4.1 


253 


1.145+0.072+0.019 


iviv 


83.8 % 


89.7% 


1137.7± 8.5 


1188 


1.044+0.030+0.007 


qqez/ 


88.3% 


93.2% 


1597.5± 9.8 


1585 


0.992+0.025+0.006 


qqjiv 


92.8 % 


96.8% 


1616. 7± 5.1 


1581 


0.978+0.025+0.003 


qqrv 


70.1% 


84.1% 


1407.8± 23.6 


1406 


0.999+0.027+0.017 


qqtv 


83.8 % 


91.7% 


4622.0± 27.6 


4572 


0.989+0.015+0.006 


qqqq 


85.9 % 


77.4% 


5845.2± 67.5 


5933 


1.015+0.013+0.012 


Total 


85.2 % 


84.7% 


11604.8± 73.4 


11693 


1.008+0.009+0.006 



Table 8: Selected events in the each of the 10 W + W~ decay topologies compared to the SM 
expectation. Also listed are the combined numbers for the six iviv decay channels and for the 
three qqlv decay channels. The efficiencies and purities for the iviv (qq^z/) decay channels 
are calculated treating all iviv (qq^z/) events as signal; e.g. the quoted efficiencies in the 
iviv channels represent the selected CC03 cross section for any iviv flavour divided by the 
generated CC03 cross section in the specific channel. Note that the total ratio of data to MC 
is for the sum of signal and background events. 



(v^)/GeV 


ff ww [pb] 


ct sm fpbl 

WW J 


161.30 


3.56+ 0.88+ 0.11 


3.61 


172.11 


12.14+ 


1.34+ 


0.22 


12.10 


182.68 


15.38+ 


0.61 + 


0.13 


15.37 


188.63 


16.22+ 


0.35+ 


0.11 


16.26 


191.61 


15.87+ 


0.86+ 


0.10 


16.55 


195.54 


18.21 + 


0.57+ 


0.12 


16.82 


199.54 


16.23+ 


0.54+ 


0.11 


17.00 


201.65 


17.94+ 


0.81 + 


0.11 


17.05 


204.88 


15.99+ 


0.52+ 


0.11 


17.10 


206.56 


17.58+ 


0.42+ 


0.12 


17.12 



Table 9: Measured CC03 W + W cross sections from a combined fit to all data. The last 
column shows the SM expectations which are taken from the average of the predictions from 
YfsWW and RacoonWW. 
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by the SM expectation scaled by a single data/SM ratio gives: 

data/SM = 1.002 ± 0.011(stat.) ± 0.007(syst.) ± 0.005(theory), 
where the SM expectation is the mean of the cross sections predicted by YfsWW and RacoonWW. 

5 Measurement of the W Branching Fractions 

A simultaneous fit to the numbers of W + W~ candidate events in the ten identified final states 
(ez/ e ez/ e , iiv^v^ tv t tv t , ei/ e /^, eu e ru T , iiv l jv T , qqeu e , qq^v^ qqru T , and qqqq) observed 
by OPAL at each of the ten centre-of-mass energies between 161 GeV and 207 GeV gives the 
following values for the leptonic branching fractions of the W boson: 

Br(W^ez/ e ) = 10.71 ± 0.25(stat.) ± 0.11(syst.) % 
Br(W^/i^) = 10.78 ± 0.24(stat.) ± 0.10(syst.) % 
Br(W^rz/ T ) = 11.14 ± 0.31(stat.) ± 0.17(syst.) %. 

Correlations between the systematic uncertainties at the different energy points have been ac- 
counted for in the fit as have correlations in the selection efficiency uncertainties for the different 
channels. These results are consistent with the hypothesis of lepton universality, and agree well 
with the SM prediction of 10.83% [5]. The correlation coefficient for the resulting values of 
Br(W — > eu e ) and Br(W — > fiu^) is +0.14. The correlation coefficients for Br(W — > ez/ e ) and 
Br(W — > fiu^) with the measurement of Br(W — > rv r ) are —0.30 and —0.23 respectively. A 
simultaneous fit assuming lepton universality gives 

Br(W^qq) = 67.41 ± 0.37(stat.) ± 0.23(syst.) %, 

which is consistent with the SM expectation of 67.51 %. Here, the largest single source of 
systematic uncertainty is that from the e + e~ — > qq background in the W + W~— > qqqq channel. 

Assuming the quark-lepton universality of the strength of the charged current weak inter- 
action, the hadronic branching fraction can be interpreted as a measurement of the sum of 
the squares of the six elements of the CKM mixing matrix, \Vij\, which do not involve the top 
quark: 

Br(W -> qq) / a s (M w ) 

(l-Br(W^qq)) \ n 

The theoretical uncertainty of this improved Born approximation due to missing higher order 
corrections is estimated to be 0.1% [5]. Taking q s (Mw) to be 0.119 ±0.002 [44], the branching 
fraction Br(W — > qq) from the 161 — 209 GeV data yields 

\Vij\ 2 = 1.993 ± 0.033(stat.) ±0.023(syst.), 

«=u,c; j=d,s,b 

which is consistent with the value of 2 expected from unitarity in a three-generation CKM 
matrix. If one assumes unitarity and a three-generation CKM matrix then this measurement 
can be interpreted as a test of quark-lepton universality of the weak coupling constant for 
quarks, g^, and for leptons, g^: 

g$/g% = 0.996 ± 0.017(stat.) ± 0.011(syst.). 
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Finally, using the experimental measurements of the CKM matrix elements other than 
\V CS \ gives |Kd| 2 + |Ks| 2 + |Kb| 2 + |Kd| 2 + |Kb| 2 = 1-054 ± 0.005 [44], and the OPAL result 
for J2i=u c -j=d s b \Vij\ 2 can be interpreted as a measurement of |V^| which is the least well 
determined of these matrix elements: 

\V CS \ = 0.969 ± 0.017(stat.) ±0.012(syst.). 

The uncertainty in the sum of the other five CKM matrix elements, which is dominated by the 
uncertainty on \V c d\, contributes a negligible uncertainty of 0.003 to this determination of |V^ S |. 

6 e + e _ — ► W + W~ Differential Cross Section 

In qqiu events it is possible to reconstruct the polar angle of the produced W~ with respect 
to the e~ beam direction, cos#w-, where the charge of the lepton tags the W ± and the jet 
momenta and the remaining event properties give the direction. Selected qqev e and qq/xz^ 
events are used to measure the differential cross section, d(t7 ww )/d(cos0 w -). Events selected 
solely by the trackless selections are not used here. Selected qqru T events are not considered 
due to the larger background and less reliable determination of lepton charge resulting from 
the possibility of the candidate tau being formed from tracks from the fragmentation of the 
quarks. 

The measured qqez/ e and qq/if^ differential cross sections are corrected to correspond to the 
CC03 set of diagrams but with the additional constraint that, at generator level, the charged 
lepton is more than 20° away from the e + e~ beam direction, 20° < 0£± < 160°. This angular 
requirement is closely matched to the experimental acceptance. It also greatly reduces the 
difference between the full four-fermion cross section and the CC03 cross section by reducing the 
contribution of t-channel single-W diagram in the qqeu e final state. At the MC generator level 
the angle cos 6> w - is defined in terms of the four-momenta of the fermions from the W~ decay 
using the CAL05 photon recombination scheme [26]. The quoted differential cross sections 
correspond to d[<r(e + e~ — > W + W~ — > qqez/ e ) + <r(e + e~ — > W + W~ — > qq/xz/ /i )]/dcos6 l w - within 
the above generator level acceptance. 

The differential cross section is measured in ten bins of cos W - with the data divided into 
four v/i ranges: 180.0-185.0 GeV; 185.0-194.0 GeV; 194.0-202.5 GeV; and 202.5-209.0 GeV. 
Experimentally the angle cos#\y- can be obtained from the measured momenta of the two jets 
with the lepton used to tag the charge of the W boson. However, to improve the angular 
resolution a kinematic fit to the four momenta of the two jets and the lepton is employed [31]. 
If the fit converges with a fit probability of > 0.1 % [31] the fitted jet momenta are used. If 
the kinematic fit yields a fit probability of < 0.1 %, which is the case for approximately 4% of 
qqtv events, cos 9 W - is calculated from the measured jet four-momenta. From MC the cos 6\y- 
resolution is found to be approximately 0.05. 

The reconstructed cos 6* w - distributions are corrected to the signal definition using the MC 
background estimates and a simple bin-by-bin efficiency correction. It has been verified that 
this simple bin-by-bin correction method is in good agreement with a more complete unfolding 
using the reconstructed to generator level migration. 

The systematic uncertainties on the selection efficiencies and background cross sections 
described above are propagated to the differential cross section measurement. In addition it 
is known from studies of lepton pair production at LEP1 that the OPAL MC underestimates 
the fraction of events where the lepton track is assigned the wrong charge [45]. This arises 
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from imperfect tracking in the region of the jet chamber anode planes. For the data considered 
here the MC predicts that 0.5 % of tracks are assigned the wrong charge. Based on previous 
studies [45] it is estimated that the corresponding number for data is (1.0 ±0.5) %. In deriving 
the efficiency corrections, the MC reconstructed cos^w- distributions are corrected for this 
difference and the full size of the correction is taken as the charge identification systematic 
uncertainty. 

The measured differential cross sections in the 10 bins of cos 9^- for the four energy ranges 
are shown in Figure 5 and the results are given in Table 10. The data are in good agreement with 
the SM expected generator level distributions obtained from either YfsWW or RacoonWW. 
Although the differential cross sections for these data have not been published previously it 
should be noted that a deviation from the SM would have shown up in the OPAL triple gauge 
coupling analysis [46] which uses similar distributions. 



bin 



(y/s) = 182.7 GeV 



Differential cross section [pb] 
(y/s) = 189.0 GeV (y/s) = 198.4 GeV 



(y/s) = 205.9 GeV 



-0.8 
-0.6 
-0.4 
-0.2 
0.0 
+0.2 
+0.4 
+0.6 
+0.8 
+1.0 



0.44 + 0.22 
0.90 + 0.30 
1.09 + 0.31 
1.24 + 0.33 
1.91 + 0.41 
2.29 + 0.45 
2.40 ± 0.46 
2.88 + 0.51 
3.87 + 0.60 
4.77 + 0.69 



±0.02 
±0.02 
±0.01 
±0.01 
±0.01 
±0.01 
±0.01 
±0.02 
±0.02 
±0.03 



0.60 ± 
0.97 ± 
1.00 ± 
1.12 ± 

1.19 ± 
1.95 ± 

2.20 ± 
2.71 ± 
3.64 ± 
5.83 ± 



0.14 ±0.03 
0.16 ±0.02 
0.16 ±0.01 
0.17±0.01 
0.17±0.01 
0.21 ±0.01 
0.23 ±0.01 
0.26 ±0.01 
0.31 ±0.02 
0.40 ± 0.04 



0.62 ± 
0.66 ± 
0.83 ± 
1.39 ± 
1.52 ± 
1.95 ± 
1.85 ± 
2.41 ± 
4.19 ± 
6.98 ± 



0.15 ±0.04 
0.15 ±0.02 
0.15 ±0.01 
0.19 ±0.01 
0.20 ±0.01 
0.22 ±0.01 
0.22 ±0.01 
0.25 ±0.01 
0.34 ±0.03 
0.47 ±0.04 



0.46 ± 
0.59 ± 
0.44 ± 
0.98 ± 
1.14 ± 
1.96 ± 
2.31 ± 
2.91 ± 
4.59 ± 
7.23 ± 



0.12 ±0.04 
0.13 ±0.02 
0.11 ±0.02 
0.15 ±0.01 
0.16 ±0.01 
0.21 ±0.01 
0.23 ±0.01 
0.26 ±0.02 
0.33 ±0.03 
0.44 ±0.05 



Table 10: The measured differential cross section, d[o"(e + e~ — > W + W~ — > qqez/ e ) + <r(e + e~ — > 
W + W~ — > qq/iz/^)]/dcos#w- expressed in ten bins of cos#\y- f° r the four centre-of-mass en- 
ergy ranges. The cross sections correspond to the CC03 set of diagrams with the additional 
requirement that the charged lepton is more than 20° from the beam axis, 20° < 9 f ± < 160°. 
For each entry, the first uncertainty is statistical and the second systematic. 



7 Conclusions 

From a total data sample of 701.1 pb -1 recorded with e + e~ centre-of-mass energies of y/s = 
161-209 GeV with the OPAL detector at LEP 11693 W-pair candidate events are selected. The 
combined data samples is almost a factor three larger than the previous OPAL publication. 
This large sample of events has enabled a significant reduction in a number of systematic 
uncertainties compared with our previous publications. 

The data are used to test the SM description of W + W~ production in the centre-of-mass 
range y/s = 161 — 209 GeV. The W-pair production cross sections at 10 different centre-of-mass 
energies are found to be consistent with the Standard Model expectation: 

data/SM = 1.002 ± 0.011(stat.) ± 0.007(syst.) ± 0.005(theory). 
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The data are then used to determine the W boson leptonic branching fractions: 



Br(W -> ev e ) 
Br(W -> /^) 
Br(W -> ri/ T ) 



10.71 ± 0.25(stat.) ± 0.11(syst.) % 
10.78 ± 0.24(stat.) ± 0.10(syst.) % 
11.14 ± 0.31(stat.) ± 0.17(syst.) %. 



These results are consistent with lepton universality of the charged current weak interaction 
and with the results of the other LEP collaborations [47-49]. Assuming lepton universality, the 
branching ratio to hadrons is determined to be 67.41 ± 0.37(stat.) ± 0.23(syst.) % from which 
the CKM matrix element \V CS \ is determined to be 0.969 ± 0.017(stat.) ± 0.012(syst.). The 
differential cross section as a function of the W~ production angle is measured for the qqez/ 
and qq/iz/ final states and found to be consistent with the SM expectation. 
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Figure 1: Distributions of (a) the total visible energy in the event scaled to the centre-of- 
mass energy, (b) the magnitude of the net visible transverse momentum in the event scaled to 
the beam energy, (c) the reconstructed total visible invariant mass of the event, and (d) the 
invariant mass of the system recoiling against the visible system. All plots show the selected 
Ivlv events for the combined sample from data recorded at \fs = 189 — 209 GeV. In (d) the 
events in the first bin are where the reconstructed recoil mass squared is negative. The data 
are shown as the points with error bars (statistical errors only). The total Standard Model MC 
prediction is shown by the unshaded histogram. The background components are also shown: 
interfering Ivlv (singly-hatched), non-interfering tvtv (cross-hatched) and two fermion/multi- 
peripheral (densely cross-hatched). The MC is normalised to the integrated luminosity of the 
data. 
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OPAL Vs=183-209 GeV 




Figure 2: Distributions of measured energies of the electrons, muons and visible tau decay 
products for events selected as qqez/, qqptv, and qqrz/ respectively. The combined distribution 
for all events selected as qqiu is also shown. The data are shown as the points with statistical 
error bars, while the histogram is the total MC expectation. The combined background from 
two-fermion and two-photon processes is shown by the cross-hatched region, while the non- 
CC03 four-fermion background is shown by the single-hatched region. 
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OPAL Vs=183-209 GeV 
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Figure 3: Distributions of the variables (described in the text) used in the likelihood selection 
of W + W~^ qqqq events (a)-(d) and the resulting relative likelihood distribution (e). All plots 
are shown for the combined sample from data recorded between y/s = 183 — 209 GeV. The data 
are shown as the points with error bars (statistical errors only). The total Standard Model MC 
prediction is shown by the unshaded histogram. The background components are also shown: 
four-fermion background (singly-hatched) and two-fermion background (cross-hatched). The 
MC is normalised to the integrated luminosity of the data. 
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Figure 4: The measured WW cross sections from fits assuming SM W decay branching fractions. 
The measured cross sections (points) are compared to the SM expectation (line) which is the 
average of the predictions from YfsWW and RacoonWW. The shaded region shows the 0.5 % 
theoretical error. 
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OPAL Vs=183-209 GeV 
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Figure 5: The measured W~ polar angle differential cross section for qqez/ and qq/jv events 
within the acceptance defined in the text. The measurements are shown for the four energy 
bins described in the text. The measured cross sections (points) are compared to the theoretical 
expectations (histograms) from YfsWW and RacoonWW (indistinguishable on this scale). 
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